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Abstract:  We  use  the  coupling  matrix  formalism  to  investigate  continuous- 
wave  and  pulse  propagation  through  microring  coupled-resonator  optical 
waveguides  (CROWs).  The  dispersion  relation  agrees  with  that  derived 
using  the  tight-binding  model  in  the  limit  of  weak  inter-resonator  cou¬ 
pling.  We  obtain  an  analytical  expression  for  pulse  propagation  through 
a  semi-infinite  CROW  in  the  case  of  weak  coupling  which  fully  accounts 
for  the  nonlinear  dispersive  characteristics.  We  also  show  that  intensity 
of  a  pulse  in  a  CROW  is  enhanced  by  a  factor  inversely  proportional  to 
the  inter-resonator  coupling.  In  finite  CROWs,  anomalous  dispersions 
allows  for  a  pulse  to  propagate  with  a  negative  group  velocity  such  that 
the  output  pulse  appears  to  emerge  before  the  input  as  in  “superluminal” 
propagation.  The  matrix  formalism  is  a  powerful  approach  for  microring 
CROWs  since  it  can  be  applied  to  structures  and  geometries  for  which 
analyses  with  the  commonly  used  tight-binding  approach  are  not  applicable. 
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1.  Introduction 


Coupled  optical  resonators  are  becoming  important  in  nonlinear  optics  research  as  well  as  in 
telecommunication  applications  [1,  2,  3,  4,  5,  6],  Systems  consisting  of  a  few  coupled  res¬ 
onators,  say  1  <  N  <  5,  have  been  proposed  for  optical  filtering  and  modulation  [7,  8,  9].  On 
the  other  extreme,  “large”  systems,  say  N  >  10,  can  be  regarded  as  a  new  type  of  waveguide 
termed  Coupled-Resonator  Optical  Waveguide  (CROW)  with  unique  and  controllable  disper¬ 
sion  properties  [3,  1,  2,  10], 

The  “large”  chains  (CROWs)  have  been  previously  analyzed  using  a  tight-binding  formalism 
[1].  In  the  tight-binding  method,  we  approximate  the  electric  field  of  an  eigenmode  E/j  of  the 
CROW  as  a  Bloch  wave  superposition  of  the  individual  resonator  modes  Eo  [1], 

Ejj(r,f)  =  £oexP('t%f)LexP(~  inKA)En(r  —  «Az),  (1) 

n 


where  the  «th  resonator  in  the  chain  is  centered  at  z  =  nA.  Under  the  assumption  of  symmetric 
nearest  neighbor  coupling,  the  dispersion  relation  of  the  CROW  is  [1] 


(Ok  =  £2 


Aoc 

1 - —  +  K\  cos(WA) 


(2) 
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where  Q.  is  the  resonant  frequency  of  an  individual  resonator  and  Aa  and  K\  are  defined  as 

A  a  =  J  d3r[e(r)-e0(r)]En(r)'En(r)  (3a) 

K\  =  J d3r[eo(r  —  Az)  —  e(r  —  Az)]  x  En(r)  -Ejj(r  —  Az).  (3b) 

e(r)  is  the  dielectric  coefficient  of  the  CROW  and  £o(r)  is  the  dielectric  coefficient  of  an  indi¬ 
vidual  resonator.  Therefore,  the  coupling  parameter  K\  represents  the  overlap  of  the  modes  of 
two  neighboring  resonators  and  Aa/2  gives  the  fractional  self  frequency  shift  of  C0k- 

Even  though  much  of  the  theoretical  work  on  CROWs  is  based  on  the  tight-binding  method 
[2,  5],  the  formalism  is  not  convenient  for  practical,  physical  systems.  For  example,  it  does 
not  account  for  input/output  coupling,  loss,  different  resonator  sizes,  or  variations  in  coupling 
strengths.  With  the  aim  of  rigorously  analyzing  realistic  CROW  structures,  we  use  a  matrix 
approach  [11,  12,  13]  to  study  a  system  consisting  of  N  coupled  ring  resonators  with  input  and 
output  waveguides.  Since  the  modal  properties  of  ring  resonators  can  be  easily  tailored  and 
their  fabrication  technology  is  mature  [14,  15],  they  may  enable  practical  implementations  of 
CROWs. 

2.  Transfer  matrix  formalism 

We  first  consider  an  infinite  chain  of  coupled  ring  resonators  in  order  to  obtain  its  dispersion 
relation.  Both  forward  and  backward  propagating  waves  exist  in  an  individual  resonator,  as 
shown  in  Fig.  1.  We  assume  the  coupling  region  is  sufficiently  long  compared  to  X,  so  that 


Fig.  1.  An  infinitely  long  chain  of  coupled  ring  resonators,  with  the  forward  and  backward 
propagating  field  components  labelled. 


the  light  circulating  in  one  direction  in  a  resonator  is  phase-matched  to  only  one  of  the  two 
degenerate  counter-propagating  modes  of  the  adjacent  resonator.  Using  the  notation  of  Fig.  1, 
the  coupling  between  two  adjacent  rings  can  thus  be  described  by  [16] 


'  K  ' 

t  K 

a'n 

■  d'n ' 

t  K 

'  c'r,  ' 

bn+l 

—  K*  t* 

&n+\ 

dn+\ 

—  K*  t* 

Cn+ 1 

where  t  and  K  are  respectively  the  dimensionless  transmission  and  coupling  coefficients  over 
the  coupling  length.  The  matrix  is  unitary  and  unimodular  so  that  \t |2  +  |  k\2  =  1.  Defining  a 
vector  with  the  different  field  components, 


a 


d 


n 


(5) 
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Equation  (4)  can  be  rewritten  as 


%n-\- 1  — 


P  0 
0  P 


=  Pxl 


1 

P=- 

K 


-t  1 
-1  t* 


(6a) 


(6b) 


As  the  field  propagates  around  the  ring,  it  accumulates  a  phase  shift  and  may  be  attenuated. 


Q  = 


0  Q 
Q  o 

o 


x„  =  ( 


&ipRn  o 


(7a) 

(7b) 


In  the  above  definition,  R  is  the  ring  radius  and  p  =  n(co)co/c+ia,  where  n(co)  is  the  frequency 
dependent  effective  index  and  a  is  the  loss  (or  gain)  per  unit  length  in  the  ring.  Combining  (6) 
and  (7),  we  have 

Xn+1  =  PQx„  (8) 


Equation  (8)  is  completely  general.  The  matrices  P  and  Q  can  be  specified  at  each  frequency 
to  account  for  any  frequency  dependence  of  the  effective  index,  loss,  and  transmission  and 
coupling  coefficients. 


3.  CROW  dispersion  relation 

From  a  theoretical  point  of  view,  it  is  important  to  understand  how  the  tight-binding  and  matrix 
approaches  are  related  to  each  other.  We  shall  show  the  matrix  method  embodied  in  Eq.  (8) 
converges  to  the  tight-binding  result  in  Eq.  (2)  under  certain  approximations.  The  approach  we 
adopt  is  similar  to  the  transfer  matrix  analysis  of  a  Bragg  stack  [11]. 

The  field  in  one  resonator  of  the  CROW  as  specified  by  xn  is 


fnw/  fl»  exp[^^(7T  -  0)]  +  d„  exp[-ipR{n  -  0)]  0  <  (j)  <  71 

\  b„exp[-iPR(n  +  4>)]  +  cnexp[ipR(n+4>)]  -tt<0<  0 


(9) 


where  0  is  the  azimuthal  angle  relative  to  the  propagation  direction  in  the  counter-clockwise 
sense,  and  p  is  the  radial  co-ordinate.  For  a  mode  of  an  infinite  chain  of  ring  resonators,  the 
fields  are  periodic  at  the  lattice  constant,  A.  So  applying  Bloch’s  theorem, 


x„+i  =  exp(— iKA)xn, 


(10) 


where  K  is  the  CROW  propagation  constant.  Combining  this  requirement  with  (8)  leads  to 

Det|PQ  — exp(— iWA)l/|  =  Det|(Pg)2  —  exp(-/2/C\)C|  =  0,  (11) 


where  U  is  the  identity  matrix. 

We  assume  lossless  propagation  and  Im()f)  ^  Re  (if)  for  phase-matched  coupling.  We  recall 
that  at  the  resonant  frequency  of  an  individual  resonator,  O,  ilii(il)R/c  =  m,  where  m  is  an 
integer,  and  n[Q.)  is  the  effective  index  at  Q.  Therefore,  approximating  n(Q)  «  h(cok),  we 
solve  Eq.  (1 1)  to  obtain 

(Or 

sin(— nm)  =  ±Im(K’)cos(/fA),  (12) 

which  is  the  desired  dispersion  relation  for  a  ring  CROW.  This  relation  is  exact  in  the  sense  that 
it  involves  no  assumption  about  the  coupling  strength. 
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If  we  expand  Eq.  (12)  in  the  parameter  A(Omn/D.,  Aft)  =  (Ok  —  H,  we  obtain  to  first  order 

-^  =  1±k-2cos(^A),  (13) 

where  Kn  =  Im(fc)  /{mri).  The  two  dispersion  relations  corresponding  to  the  ‘±’  coexist  for 
an  infinite  structure  to  allow  for  both  forward  and  backward  wave  propagation  (i.e.  positive 
and  negative  group  velocities).  Physically,  for  a  finite  structure  without  reflection  and  a  uni¬ 
directional  input  as  in  Fig.  3,  only  the  dispersion  relation  with  the  matching  group  and  phase 
velocities  as  the  input  wave  will  be  of  significance. 

Equation  (1 3)  is  of  a  form  identical  to  the  tight-binding  result  in  Eq.  (2).  The  correction  Aa /2 
term  does  not  explicitly  appear  in  Eq.  (13)  since  it  is  accounted  for  by  Re(x').  From  Eq.  (13), 
it  follows  that  for  Acomn/Q.  <C  1,  it  is  necessary  that  x*  <C  1.  This  condition  and  the  absence 
of  all  but  the  nearest  neighbor  coupling  are  thus  the  validity  conditions  for  the  tight-binding 
approximate  result  (13). 


Fig.  2.  The  exact  and  cosine-approximate  (i.e.  tight-binding-approximate)  dispersion  rela¬ 
tions  for  m  =  100  and  K  =  —0.8/. 


Figure  2  shows  the  dispersion  relations  for  K  =  —0.8 i  and  m  =  100  as  calculated  using  the 
“exact”  form  in  Eq.  (12)  and  the  approximated  form  in  Eq.  (13).  As  (Ok/Q  increases,  the  exact 
dispersion  relation  deviates  more  significantly  from  the  cosine  form.  For  smaller  values  of  K, 
the  deviation  from  the  cosine  dispersion  relation  is  reduced. 

The  Bloch  modes  of  the  CROW  are  given  by  the  eigenvectors  of  PQ.  At  each  frequency,  there 
are  4  Bloch  modes  corresponding  to  the  4  eigenvalues  (i.e.  values  of  K ).  The  eigenvalues  are 
exp(-iWiA)  =  <gi,  exp[— /(Wi  A+  7r)]  =  -|i,  exp(-;W2A)  =  £2,  and  exp[-t(^T2A-|-^)]  =  -<*2. 
The  corresponding  (un-normalized)  eigenvectors  are 


= 


where 


"C  +  7" 

-(C  +  7)' 

"C-7 

'-(C-7)' 

1 

-1 

1 

-1 

C  +  7 
1 

>  <Mt  = 

C  +  7 

1 

>  ki  = 

C-7 

1 

>  = 

C-7 

1 

7  = 


2 1 


,  ,—4 imnco.  „  ,~2imKC0 W1  „  ,. 

1  +exp( - — - )  +2exp( - — - )(1  -2 12), 


?  =  2t 


,  .  2imna> , 

1  +  exp(  —  ) 


,  (14) 

(15a) 

(15b) 
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The  4  eigenvectors  are  orthogonal  to  each  other  and  they  represent  standing  waves  in  each 
resonator.  In  the  limit  of  weak  coupling  |  k\  <C  1  and  ft)  «  Q,  such  that  7  «  j  k|  «  0,  £  «  1,  and 
Ci  =  ^2.  the  4  eigenvectors  reduce  to  2  degenerate  eigenvectors,  representing  the  two  different 
superpositions  of  the  clockwise  and  counter-clockwise  propagating  waves  in  a  single  resonator: 


'l' 

'-1' 

1 

-1 

kl  =  kl  = 

1 

>  = 

1 

1 

1 

In  the  limit  of  strong  coupling  and  fi)«il,t<l,  7+  C,  ss  ~  7,  and  £  —  7 

The  eigenvectors  become 


■T 

t 

" _ 2" 

t 

'o' 

'  0  ' 

1 

-1 

i 

-1 

k  = 

2 

>  = 

2 

>  ki  = 

0 

-  = 

0 

1 

1 

1 

1 

0. 


(17) 


We  observe  that  2  field  components  are  significantly  stronger  than  the  other.  This  corresponds 
to  a  wave  that  “zig-zags’'  through  the  resonators  without  making  complete  round-trips  in  each 
resonator.  The  asymptotic  behavior  of  the  eigenvectors  confirms  the  physical  picture  that  as 
|  K"|  — >  0,  the  modes  of  the  CROW  are  essentially  the  modes  of  the  independent  resonators,  and 
as  |  K"|  — s-  1,  the  microrings  no  longer  act  as  resonators  and  the  CROW  modes  are  essentially 
conventional  waveguide  modes. 

4.  Finite  CROWs  and  a  travelling  wave  picture 

For  physical  realizations  of  CROWs,  we  are  interested  in  finite  structures  with  input  and  output 
coupling.  These  properties  cannot  be  readily  dealt  with  using  the  tight-binding  method,  but  can 
be  easily  incorporated  into  the  transfer  matrices.  Figure  3  shows  a  typical  implementation  of 
a  microring  CROW:  light  is  coupled  into  and  out  of  a  set  of  coupled  ring  resonators  via  the 
input  and  output  waveguides.  Assuming  that  the  coupling  length  between  waveguides  and  the 
CROW  is  long  compared  to  A,  then  only  the  travelling  wave  phase-matched  to  the  input  can  be 
excited. 


Fig.  3.  A  CROW  consisting  of  N  ring  resonators  with  input  and  output  waveguides. 


Adopting  the  notation  in  Fig.  3,  the  fields  between  adjacent  resonators  are  related  by 


(18) 
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where  P  and  Q  are  defined  in  Eq.  (6)  and  Eq.  (7). 

By  cascading  the  transfer  matrices,  PQ,  we  obtain  an  expression  for  the  field  components  at 
the  output  of  the  CROW  after  N  identical  rings: 


aN+ 1 
pN+ 1 

=  PoutQ{PQ)N~1Pin 

1  1 

o  o 

i _ i 

= 

A 

C 

B' 

D 

ao 

bo 

where  Pm  and  Pout  describe  the  coupling  between  the  CROW  and  the  input/output  waveguides. 
For  a  single  input  to  the  waveguide,  we  set  a^+ 1  =  0.  Therefore,  the  transfer  functions  at  the 
“through”  and  “output”  ports  as  shown  in  Fig.  3  are 

—  =  ~=Ttkr{a),  (20a) 

do  B 

—  =C-^-  =  Tout{co).  (20b) 

do  D 

As  in  microring  filter  design  [17],  the  coupling  between  the  waveguides  and  the  CROW  can  be 
selected  to  maximize  the  flatness  of  the  transmission  response.  Therefore,  a  finite  CROW  can 
be  designed  to  mimic  an  infinite  CROW  over  a  bandwidth  with  a  sufficiently  flat  transmission 
response. 

An  advantage  of  the  matrix  formalism  is  that  it  is  valid  for  chains  of  any  length  N,  which 
is  essential  in  analyzing  any  physical  realization  of  a  CROW.  From  the  phase  response  of  the 
transmission  function  given  by  Eq.  (20),  we  can  deduce  the  dispersion  relation  of  the  structure. 
However,  we  note  that  the  travelling  wave  is  not  an  eigenmode  of  the  CROW,  since  the  Bloch 
modes  as  given  by  the  eigenvectors  of  PQ  are  standing  waves  as  in  Eq.  (14).  A  travelling  wave 
solution  is  formed  by  a  superposition,  either  the  sum  or  difference,  of  the  two  Bloch  modes 
with  equal  group  velocities  (i.e.  q ^  and  q  ^ ,  or  q^2  and  q  ).  The  travelling  wave  is  an  eigen¬ 
vector  of  {PQ)2,  and  it  is  verified  that  the  sense  of  propagation  in  the  rings  alternates  between 
clockwise  and  counter-clockwise  with  each  operation  of  PQ,  as  depicted  in  Fig.  3.  Therefore, 
taking  the  phase  difference  accumulated  over  two  rings  to  be  —IKK,  where  K  is  the  Bloch 
wave  vector,  such  that  the  phase  difference  between  the  output  and  the  input  is  approximately 
—  (A  —  1)AA,  we  can  determine  the  CROW  dispersion  from  the  finite  structure. 

As  an  example,  we  compute  the  dispersion  relation  of  a  finite  CROW  consisting  of  20  cou¬ 
pled  rings  with  inter-resonator  and  waveguide-resonator  coupling  constants  of  —0.5/.  The  rings 
are  lossless,  and  their  radius  is  16.4/xm.  ne g-  is  taken  to  be  constant  and  equal  to  1.5.  Figure  4 
compares  the  dispersion  relation  extrapolated  from  the  finite  CROW  with  the  dispersion  rela¬ 
tion  of  an  infinite  CROW  as  given  by  Eq.  (12).  The  small  amplitude  ripples  are  manifested  at 
the  resonance  frequencies  of  the  finite  structure.  In  the  limit  of  an  infinite  number  of  resonators, 
the  resonance  peaks  will  be  infinitesimally  close  to  each  other  and  the  ripples  will  be  smoothed 
out. 

5.  Pulse  propagation 

Pulse  propagation  through  CROWs  are  of  particular  technological  interest,  since  information 
transmitted  in  optical  communication  systems  is  typically  encoded  in  pulses.  Using  the  results 
from  the  previous  sections,  we  can  analytically  and  numerically  study  optical  pulse  propagation 
in  semi-infinite  and  finite  microring  CROWs. 

5.1.  Semi-infinite  case 

A  semi-infinite  microring  CROW  consists  of  an  infinitely  long  CROW  coupled  to  a  single  input 
waveguide  as  in  Fig.  (5).  The  input  waveguide  ensures  that  only  a  pulse  of  positive  (or  negative) 
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Fig.  4.  The  exact  dispersion  relation  for  an  infinite  CROW  and  the  dispersion  relation  as 
extracted  from  20  coupled  resonators.  The  rings  have  a  radius  of  16.4/im  and  the  inter¬ 
resonator  coupling  is  —0.5 i. 


group  velocity  propagates  through  the  structure.  Assuming  that  the  bandwidth  of  the  input  pulse 
is  within  the  bandwidth  of  the  CROW  band  such  that  all  of  the  input  light  is  coupled  into  the 
waveguide,  the  field  amplitude  b\  in  the  first  resonator  is  b\  (ft))  =  — 1/ Ki„ao(co),  where  K"„,  is 
the  coupling  coefficient  between  the  input  waveguide  and  the  first  resonator.  Since  |  Kin  |  <  1, 
the  intensity  of  the  field  inside  the  CROW  is  higher  than  that  of  the  input  pulse  by  l/|fc,„|2. 
This  does  not  violate  energy  conservation,  as  the  increased  intensity  is  a  consequence  of  the 
reduced  group  velocity  and  hence  the  spatial  compression  of  the  pulse  inside  the  CROW.  Using 
the  dispersion  relation  in  Eq.  (13),  the  maximum  group  velocity  in  the  CROW  is 


U?,  max  — 


\k\AQ. 


mn 


Defining  the  “slowing”  factor  as  in  Ref.  [18]  to  be 

S  = - — 


^effU’.max 

and  approximating  A  ~  2 R,  S  can  be  expressed  as 

n 


S  = 


2\k\ 


(21) 


(22) 


(23) 


Therefore,  for  1 qn  =  K,  the  intensity  inside  the  rings  is  roughly  enhanced  by  (|5')2.  This  result 
makes  intuitive  sense  since  the  only  loss  mechanism  for  the  otherwise  lossless  resonators  is 
the  inter-resonator  coupling.  Interestingly,  even  though  the  energy  velocity  of  the  Bloch  modes 
at  Q  corresponds  to  the  group  velocity  v^?max[l  1],  the  energy  velocity  of  a  wave  that  is  fully 
coupled  into  the  semi-infinite  CROW  is  proportional  to  |k'|2.  Hence,  the  intensity  enhancement 
is  proportional  to  the  energy  velocity  reduction  rather  than  the  group  velocity  reduction. 

Also,  in  contrast  to  CROWs,  for  other  coupled  resonator  structures  where  there  is  no  feed¬ 
back  between  the  resonators,  such  as  the  SCISSOR  [19],  the  slowing  factor  is  approximately 
proportional  to  1  / 1 7c| 2  in  the  case  of  weak  coupling.  However,  a  CROW  has  the  advantage  that, 
even  in  the  presence  of  loss,  it  is  most  transmitting  for  the  frequencies  of  the  CROW  band, 
while  a  side-coupled  resonator  is  most  attenuating  near  the  resonant  frequency  of  the  resonator. 

To  analyze  the  temporal  dynamics  a  pulse  launched  into  the  semi-infinite  CROW,  we  adopt 
a  method  of  analysis  that  is  analogous  to  pulse  propagation  in  conventional  waveguides  such 
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as  optical  fibers.  We  shall  find  in  the  limit  of  weak  coupling,  such  that  Eq.  (13)  is  a  good 
approximation,  there  exists  a  closed-form  solution  to  the  evolution  of  any  arbitrary  input  pulse. 

Through 


Input 


Fig.  5.  A  semi-infinite  CROW. 

The  electric  field  at  where  b\  is  taken,  S'(t.z  =  0),  can  be  expressed  as  the  Fourier  integral 

g(t,z  =  0)  =  /  dcof?i(co)exp(/cur),  (24a) 

J  band 

r  dt' 

b\((0)=  /  —  &(t',z  =  0)exp(— im').  (24b) 

J  2  n 

At  z  =  A'A,  each  frequency  component,  b\  (ft)),  acquires  a  phase  shift  of  NKA,  so  the  field  is 

r  cp' 

xp(/fi)f)  / — g(t' ,z  =  0)exp[— i(cot' +  K(co)NA)].  (25) 

J  2k 


g(t1z  =  NA)  =  f  < 

J  band 


dcoexp 


However,  K(co)  is  given  by  the  dispersion  relation  of  the  CROW,  Eq.  (13).  Therefore,  instead 
of  integrating  over  frequency  in  (25),  if  we  integrate  over  the  half  of  the  Brillouin  zone  that 
gives  the  appropriate  group  velocity  (for  example,  the  right  half),  we  obtain 


£{t,z  =  NA)  =  -A  Q.K2eint  f(^/A  dK  sm(KA)e-‘KNAe‘nK2C0^KA^-^ 


(26) 


f£g(t',z  =  0)e 


-idt' 


Equation  (26)  can  be  further  simplified  by  letting  x  =  KA,  and  invoking  the  Jacobi-Anger 
expansion  [20], 

eiO.Kcos(x)(l-t')  =^CmJm[aK2(t-tr)}  cos  (mx) 


Cm  — 


1  if  m  =  0 

2 im  if  m  >  0 


to  arrive  at 

g(t,z  =  NA)  =  —  Q.K2ein‘YtmcmJo  dxsin(x)cos(mx)e~,xN 

J^Jm[aK2(t-t'))g(t',z  =  o)e-™'. 

However,  am^  =  djc sin (jc)  cos (mx)e~lxN  ^  0  only  for  certain  values  of  m  and  N: 


(27a) 

(27b) 

(28) 


UL 
'  4 


-2(m2+N2-l) 
(m2+N2~l)2—4m2N2 


for  N  =  m—  1  and  N  =  —m  —  1 
for  N  =  m  +  1  and  N  =  —  m  + 1 
for  N  +  m  =  even 


(29) 
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So  the  equation  for  the  pulse  envelope  E(t,z),  such  that  ${t,z)  =  E(t,z)eini,  is  given  by  the 
convolution  integral 

E{t,z  =  NA)  =  -^-Tcmam:N  [ dt'Jm[£lK2(t-t')]E{t',z  =  0).  (30) 

2  n  m  J 

The  Fourier  transform  of  a  Bessel  function  Jn(t)  is  only  defined  within  \2nf\  <  1  [21],  which 
accounts  for  the  finite  bandwidth  of  the  CROW,  0(1-  K3 1 )  <  ft)  <  Of  1  +  |  Ki  | ) . 

Equation  (30)  holds  for  an  arbitrary  input  pulse  and  its  sole  assumption  is  the  cosine- 
approximate  dispersion  relation,  which  is  valid  for  small  K.  The  nonlinear  dispersive  nature 
of  the  CROW  is  embodied  in  the  summation  over  the  Bessel  functions.  Figure  6(a)  shows  the 
evolution  of  a  Gaussian  input  pulse  E(t,z  =  0)  =  exp(— t2 /T2)  as  calculated  using  Eq.  (30). 
Figure  6(b)  shows  the  numerical  results  obtained  from  the  transfer  matrices.  The  analytical  so¬ 
lution  is  in  excellent  agreement  with  the  fully  numerical  approach.  As  the  pulse  propagates, 
even  though  the  main  peak  travels  at  the  group  velocity,  the  ripples  develop  only  at  the  tail  end 
of  the  pulse. 


Fig.  6.  Evolution  of  a  2.4ps  (FWFIM)  Gaussian  pulse  centered  about  1.5/im  in  a  semi¬ 
infinite  CROW  with  K %  =  0.0016.  The  fields  are  normalized  to  the  maximum  field  ampli¬ 
tude  in  the  first  resonator,  (a)  Theoretical  results  computed  using  Eq.  (30).  (b)  Results  com¬ 
puted  numerically  with  the  transfer  matrices  using  a  chain  of  100  ring  resonators  («eff  =  1-5, 
R  =  16/rm). 


5.2.  Finite  case 

Pulse  propagation  through  finite  CROWs  can  be  easily  analyzed  using  the  transfer  matrices 
results  of  Eq.  (20).  Since  Eq.  (20)  is  specified  at  each  frequency,  we  simply  have  to  find  the 
product  between  the  transfer  functions  and  the  spectral  components  of  the  input  pulse.  The 
temporal  behavior  follows  naturally  from  the  Fourier  transform. 

Distortionless  propagation  through  an  arbitrary  finite  CROW  can  always  be  achieved  if  the 
input  pulse  is  sufficiently  narrow-band  such  that  the  transmission  function  of  the  drop  port,  as 
defined  in  Fig.  3,  over  the  pulse  bandwidth  is  near  unity.  However,  short  pulses  which  become 
distorted  as  they  propagate  in  the  CROW  are  also  of  fundamental  interest.  For  this  purpose,  we 
take  an  example  consisting  of  10  coupled  ring  resonators  of  radius  164.5/im  and  neff  =  1.5. 
The  inter-resonator  coupling  constant  is  —0.3/  and  the  coupling  between  the  waveguides  and 
CROW  is  —0.5/.  The  transfer  characteristics  of  this  structure  are  shown  in  Fig.  (7).  We  launch 
a  30.5ps  (FWHM)  long  pulse  centered  at  1.55/rm  into  the  CROW. 
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Fig.  7.  The  transmission  characteristics  of  a  10  ring  long  CROW.  The  ring  radius  is 
164.5/rm  and  ne g  =  1.5.  Inter-resonator  coupling  is  —0.3/  and  the  waveguide-CROW  cou¬ 
pling  is  —0.5/.  A  30.5  ps  (FWHM)  long  pulse  centered  at  1 ,55/im  is  input  into  the  CROW, 
(a)  Transmittance  of  at  the  drop  port.  The  dashed  line  shows  the  spectrum  of  the  input 
pulse,  (b)  Phase  response  at  the  drop  port,  (c)  Transmittance  at  the  through  port,  (d)  Phase 
response  at  the  through  port. 


Using  the  transfer  matrices,  we  can  examine  how  a  pulse  evolves  in  the  CROW  by  finding  the 
transfer  functions  associated  with  an  or  bn.  Figure  8  shows  the  evolution  of  the  pulse  through 
the  CROW.  Even  though  the  output  pulse  is  attenuated  compared  to  the  input,  the  field  intensity 
inside  the  rings  can  be  greater  than  the  input,  as  in  the  case  of  the  semi-infinite  CROW.  The 
intensity  build-up  is  verified  by  a  FDTD  simulation  discussed  in  Sect.  5.3.  The  significant 
increase  in  the  intensity  of  the  input  pulse  inside  the  CROW  can  be  used  to  enhance  the  strengths 
of  nonlinear  optical  interactions.  As  noted  earlier,  we  can  account  for  loss  (or  gain)  in  our  model 
by  including  an  imaginary  part  to  the  propagation  factor  /j .  We  have  found  the  transfer  matrices 
give  excellent  agreement  with  experimental  results  [22]. 

Another  interesting  effect  is  the  small  amplitude  ripple  that  follows  the  main  peak  in  each 
resonator.  The  ripple  is  travelling  from  the  end  of  the  CROW  back  to  the  start  at  approximately 
the  group  velocity  of  the  forward  moving  pulse.  This  is  analogous  to  a  reflection  from  the  end  of 
a  waveguide,  though  in  the  microring  CROW  described  here  there  are  no  reflection  mechanisms 
as  the  coupling  is  assumed  to  be  perfectly  phase-matched.  Indeed,  Fig.  9  shows  that  the  ripple  at 
the  through  port  is  delayed  from  the  the  drop  port  pulse  by  the  travelling  time  between  the  input 
and  the  drop.  Therefore,  although  the  microring  CROW  is  composed  of  “microscopic,”  discrete 
elements,  it  possesses  certain  “macroscopic”  properties  that  mimic  conventional  waveguides. 

At  the  through  port,  we  may  obtain  negative  group  velocities,  which  some  researchers  refer 
to  as  “superluminal  propagation.”  [23,  6]  In  the  time  domain,  the  main  (highest)  peak  of  the 
output  pulse  does  indeed  appear  before  the  peak  of  the  input  pulse.  Figure  9  shows  the  output 
pulses  at  the  through  and  drop  ports  as  well  as  in  the  input  pulse.  The  peak  of  the  through  pulse 
is  approximately  5ps  before  than  the  peak  of  the  input,  as  though  the  output  appears  before 
the  input.  However,  the  pulse  is  attenuated  and  distorted.  This  behavior  is  accounted  for  by  the 
anomalous  dispersive  properties  at  the  through  port  in  Fig.  7(d).  The  anomalous  dispersion  is 
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Fig.  8.  The  pulse  transmission  through  the  CROW  described  in  Fig.  (7).  The  Otli  resonator 
is  the  input  pulse  and  the  1  Ith  resonator  is  the  output  pulse  at  the  drop  port. 


also  confirmed  by  the  FDTD  simulation  discussed  in  the  next  section. 


Fig.  9.  The  input  pulse  and  the  output  pulses  at  the  drop  and  through  ports  of  the  the  CROW 
described  in  Fig.  (7).  The  solid  vertical  line  marks  the  maximum  of  the  input  pulse,  and  the 
dashed  vertical  line  marks  the  maximum  of  the  output  pulse  at  the  through  port.  The  peak 
of  the  through  port  pulse  occurs  about  5ps  sooner  than  the  peak  of  the  input. 


5.3.  FDTD  simulations 

As  a  test  for  the  transfer  matrix  method  and  a  confirmation  of  the  intensity  build-up  and  anoma¬ 
lous  dispersion,  we  use  a  finite  difference  time  domain  (FDTD)  simulation  to  study  the  pulse 
propagation  through  two  coupled  ring  resonators.  The  waveguides  and  rings  are  0.2 /.mi  wide. 
They  are  set  in  air  and  have  an  index  of  refraction  of  3.5.  The  rings  have  a  radius  of  5/im, 
and  the  wavelength  dependent  effective  index,  as  extrapolated  from  a  separate  FDTD  simula¬ 
tion  of  the  waveguides,  is  neff  =  3.617  —  0.5539A.  The  coupling  between  the  rings  is  —0.32 i, 
and  the  coupling  between  the  rings  and  the  waveguides  is  —0.4/.  A  2.4ps  (FWHM)  Gaussian 
pulse  is  launched  into  the  system,  and  the  fields  at  the  through  port,  drop  port,  and  inside  the 
rings  are  monitored.  We  compare  the  transfer  matrix  method  with  the  FDTD  simulation  in  Fig. 
10(a),  showing  that  the  approaches  are  in  excellent  agreement.  The  anomalous  dispersion  at 
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the  through  port  and  the  increase  in  intensity  in  the  coupled  rings  are  confirmed  by  the  FDTD 
simulation  and  are  evident  in  Fig.  10(b). 


Fig.  10.  FDTD  simulation  of  2  coupled  ring  resonators  with  input  and  output  waveguides. 
The  radius  of  the  rings  is  5flm,  and  the  effective  index  is  ne ff  =  3.617  —  0.5539A.  The 
inter-resonator  coupling  is  —0.32 i  and  the  waveguide-resonator  coupling  is  —0.41.  The 
input  pulse  is  a  2.4ps  (FWHM)  Gaussian  centered  at  1.55/im.  (a)  Comparison  between 
the  FDTD  simulation  and  the  transfer  matrix  method.  Output  refers  to  the  drop  port,  (b) 
Intensity  built-up  and  anomalous  dispersion  as  confirmed  by  the  FDTD  simulation. 


6.  Comparison  with  Fabry-Perot  resonators 

The  transfer  matrix  analysis  discussed  here  is  general  and  can  be  applied  to  Bragg  stacks  and 
Fabry-Perot  resonators  as  well.  A  Bragg  stack  can  be  regarded  as  a  chain  of  coupled  Fabry- 
Perot  etalons.  Bayindir  et  al.  have  recently  demonstrated  a  CROW  operating  near  600nm  using 
coupled  Fabry-Perot  resonators  fabricated  in  a  multi-layer  fashion  [24].  There  have  also  been 
extensive  studies  in  the  linear  and  nonlinear  optical  propagation  in  Bragg  stacks  [11,  25,  26,  27, 
28,  29].  In  Fabry-Perot  resonators,  the  coupling  parameter  K  corresponds  to  the  transmission 
coefficient,  while  t  corresponds  to  the  reflection  coefficient.  They  can  be  calculated  from  the 
Fresnel  coefficients  or  from  an  analysis  of  a  Bragg  stack  in  the  case  of  a  Fabry-Perot  with 
Bragg  end  mirrors  [11].  The  two  major  differences  between  the  Fabry-Perot  resonator  and  the 
ring  resonator  are  1)  the  former  can  be  fully  described  by  a  2  x  2  transfer  matrix  while  the 
latter  requires  a  4  x  4  matrix,  and  2)  the  reflection  and  transmission  coefficients  may  contain 
a  real  part.  Therefore,  in  contrast  to  the  icos(k'A)  dependence  in  Eq.  (13),  there  is  only  one 
dispersion  curve  for  a  Fabry-Perot  CROW  at  a  given  frequency  range  and  there  may  also  be  an 
additional  phase  shift  in  the  cosine  dependence. 

In  contrast  to  ring  resonators,  the  coupling  between  Fabry-Perot  resonators  with  Bragg  end 
mirrors  is  controlled  by  Bragg  reflection.  Even  though  the  coupling  coefficient  may  be  more 
stringently  controlled  in  these  structures,  ring  resonators  remain  an  attractive  option  for  CROWs 
in  planar  integrated  optical  circuits  because  they  can  be  fabricated  in  a  single  lithographic  step. 
Moreover,  recent  developments  in  coupled  ring  resonators  in  polymer  and  semiconductor  mate¬ 
rials  illustrate  the  potential  of  using  ring  resonators  as  constituent  elements  in  CROWs  [22,  30]. 
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7.  Conclusion 


In  summary,  the  transfer  matrix  method  is  used  to  analyze  microring  coupled  resonators.  The 
transfer  matrix  and  tight-binding  approaches  yield  equivalent  dispersion  relations  in  the  limit 
of  weak  coupling.  We  also  study  pulse  propagation  through  semi-infinite  and  finite  CROWs  to 
find  intensity  enhancement  as  well  as  anomalous  dispersion.  The  matrix  method  can  account 
for  finite  chains,  holds  for  any  coupling  strength,  applies  to  travelling  waves,  and  can  treat 
heterogeneous  chains  consisting  of  an  arbitrary  mix  of  resonators  and  coupling  constants.  These 
features  make  the  transfer  matrices  versatile  for  device  design  and  for  analyzing  experimental 
results  of  microring  CROWs  [22]. 
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